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Helices are a key folding motif in protein structure. The question which factors deter- 
mine helix stability for a given polypeptide or protein is an ongoing challenge. Here we use 
van der Waals corrected density-functional theory to address a part of this question in a 
bottom-up approach. We show how intrinsic helical structure is stabilized with length and 
temperature for a series of experimentally well studied unsolvated alanine based polypep- 
tides, Ac-Alan-LysH+. By exploring extensively the conformational space of these molecules, 
we find that helices emerge as the preferred structure in the length range n=4-8 not just 
due to enthalpic factors (hydrogen bonds and their cooperativity, van der Waals disper- 
sion interactions, electrostatics), but importantly also by a vibrational entropic stabilization 
over competing conformers at room temperature. The stabilization is shown to be due to 
softer low- frequency vibrational modes in helical conformers than in more compact ones. 
This observation is corroborated by including anharmonic effects explicitly through ah initio 
molecular dynamics, and generalized by testing different terminations and considering larger 
helical peptide models. 



I. INTRODUCTION 



Polypeptide helices are a key secondary structure motif in a wide range of proteins [T-^. It is 
well known that some amino acids (e.g., alanine) exhibit a stronger helix propensity than others 
[IHH], but the fact that the helical structure is so abundant in proteins is still intriguing. From a 
thermodynamic point of view, there are at least two possible limits in which helices compete with 
other structure prototypes. Towards high temperature, one expects the transition to a random coil 
[T2] . which should become entropically favored as the temperature increases. Towards low tem- 
perature, however, helices may themselves compete with other enthalpically stable conformations. 
However, in the most interesting regime, namely intermediate, physiological temperatures, stability 
may be determined by a delicate balance between enthalpy and entropy [13j. We here unravel this 
balance quantitatively for the emergence of helical structure in a particularly well studied series 
of unsolvated polyalanine based peptides Ac-Alari-LysH+, n=4-8 [HlEJ- We consider exphcitly 
not just the helical part of conformational space, but actually the much larger, general low-energy 
conformational space of the peptides, of which helices are a part. In this paper we show: (i) a 
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comprehensive search of the conformational space for Ac-Alar2,-LysH+, n=4-8; (ii) harmonic free 
energy calculations for several structural candidates; (iii) outlook on the role of anharmonicities 
in the potential energy surface; and (iv) a theoretical comparison for longer model peptides, con- 
sidering only helical motifs, with a different termination, in order to clarify the impact of Lys 
on the soft vibrational modes. Our key finding is that there is a significant vibrational entropic 
stabilization of helices compared to other, more compact conformers, a contribution that should 
indeed make a difference in actual proteins as well. This contribution is intrinsic to the helix and 
should therefore act largely independently, not entangled with environment-dependent terms such 
as a solvent entropy. Interestingly, an essential role of low-frequency modes is also actively debated 
in other areas of protein science, related to their function [T61 - I20] . 

Beginning with the terms that shape the potential-energy surface (PES), known reasons for helix 
stability include [21, 22j (i) their efficient hydrogen-bond (H-bond) network and increasing H- 
bond cooperativity with helix length [23-.-.-27J, (ii) suitably bonded and/or electrostatically favorable 
termination, for instance the LysH^ termination considered here \14\ [28H33] , and (iii) remarkably, 
a rather specific favorable contribution of van der Waals (vdW) interactions for a-helices [271 ES] . 
Clearly, the peptide chain length plays a role: Too short chains have too few and too weak hydrogen 
bonds for helices to compensate the cost of strain in the backbone ^5] [321 [36H38] . In practice, 
environment effects will necessarily influence helix stability [ 38l441j . In an aqueous medium, the 
hydrophobic effect will be prominent [42j. In water-poor conditions, like membranes, helices are 
frequently observed [l3l W4:\. Finally, in vacuum conditions, the longer members (n >8) of the 
polyalanine based peptide series studied here assume helical structure in experiment [m HHl 26] . 
These helices are stable in vacuo even up to extreme temperatures (not expected in solution) 
[271 147] , or after soft landing on a surface |48j . 

The potential energy surface also shapes entropy, and thus the effect of temperature T. With in- 
creasing T, the conformational entropy of the backbone will favor an unfolded state [l2 l l22 l H9 l [50] 
(so-called "random coil"), while at low T helices may also compete with other, enthalpically more 
stable conformers. For instance, gas-phase ion mobility spectrometry (IMS) by Jarrold and cowork- 
ers [13j showed that the Ac-Ala4-Gly7-Ala4H+ polypeptide is helical at r=400 K but globular at 
room temperature. A similar structural change was observed in experiments involving multiply 
protonated polyalanine in the gas-phase in Ref. [51j. Empirical force- field based simulations by Ma 
and coworkers [52j of more than 60 small peptides indicate that the vibrational entropy (harmonic 
approximation) could stabilize a-helices or /3-hairpins over competing low-temperature conformers. 
Very recently, Plowright and coworkers |53| used density-functional theory (DFT) including dis- 
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persion contributions (the B97-D [54] exchange-correlation functional) to suggest that, for a small 
neutral four-residue peptide, /3-sheets and conformers containing 3io helical loops are stabilized by 
the harmonic vibrational entropy at finite temperatures. 

In the present work, we provide independent, unambiguous, and quantitative computational evi- 
dence that the vibrational entropy acts to stabilize helical conformers with increasing temperature 
over more compact, enthalpically competitive structures. The reason, in short, is traced to the 
softer low- frequency modes of helices, which are also reflected in the dynamics (anharmonic case). 
We focus on polyalanine-based peptides, since alanine is known to have a high helix propensity 
both in solution [2TJ [55] and in vacuo [8J. For Ac-Ala^-LysH+ (n=4-20) in the gas phase, IMS 
[T5] and first-principles calculations compared to experimental vibrational spectroscopy at room 
temperature jjG] suggest a cross-over from non-helical to helical preferred conformers as a function 
of polyalanine chain length. For n=5 there is a competition between different conformers, while 
the n=10 and n=15 conformers are found to be firmly in the helical range [46j. In a previous pub- 
lication [27j we have quantified, from first principles, the contributions from electrostatics, H-bond 
cooperativity, and van der Waals interactions on the stability of unsolvated polylalanine-based 
helices against unfolding. This class of systems is thus an ideal testing ground to clarify the struc- 
tural competition of non-helical (compact) and helical conformers as a function of chain length 
also toward the opposite temperature limit (low temperature, folded state). In the following, we 
address conformational preference of Ac-Ala^-LysH+ in vacuo for n=4-8, i.e., the length range in 
which the helical preference at room temperature develops. 

Our work is based on an exhaustive prediction of low-energy conformers using DFT and the PBE 
[56] exchange-correlation functional, corrected to account for long-range vdW interactions [57j (here 
called PBE+vdW). This level of theory treats accurately, and without system-specific empirical 
parameters, critical length-dependent contributions such as H-bond cooperativity [23l [251427] and 
vdW interactions [27], including their effect on vibrational frequencies and in ab initio molecular 
dynamics. 

II. METHODS 

Our conformational search strategy, used to find structure candidates of Ac-Alar^LysH^, n=4-8, 
consists of two steps. For a detailed discussion of the search strategy we use here, we refer the 
reader to Ref. [58j. Both the Lys residue and the C-terminal COOH are considered protonated 
throughout, as is known to be the gas-phase preference when the N-terminus is capped. [141 l59l [60] 
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In the first step, we begin by an extensive, unconstrained force field (OPLS-AA [61]) basin- hopping 
search (using the TINKER package [62]), aiming to simply "list" as many different conformers as 
possible. The same strategy was employed in Ref. [16], enumerating a huge number of candidate 
structures: at least 10^ conformers for each of the molecules in question. Our particular choice of the 
OPLS-AA force field is not motivated by any other reason than that an input structure "generator" 
for DFT was needed. In particular, care was taken that small changes to the parameters in the 
OPLS-AA part of the search do not affect the structure manifold considered in the DFT step 
below. [58] 

In the second step, a wide range of conformers suggested by the force field is fully relaxed using 
DFT with the PBE+vdW exchange-correlation functional in all-electron total energy calculations 
(FHI-aims program package [63]). We employ essentially converged numeric atom-centered basis 
sets and other numerical settings [63l [64] for the relaxations. In detail, the PBE+vdW part of 
our searches covers 1068, 1000, 800, 800, and 820 conformers for n=4, 5, 6, 7, 8, respectively. For 
n=6, 7, and 8, these numbers are chosen to ensure that at least all conformers identified in the 
lowest 7 kcal/mol 0.3 eV) of the force field step are included in the PBE+vdW relaxations. For 
n=4 and 5, we explored the limits of our search strategy [58j. When comparing locally relaxed 
minima in PBE+vdW that were obtained by starting from an OPLS-AA local PES minimum 
structure, we observe a correlation of the relative energy hierarchies, but with a large scatter (max. 
around 50 meV per residue). Therefore, a large number of conformers must be considered for post- 
relaxation in PBE+vdW. The comparison also reveals structure-specific force field errors that might 
otherwise go unnoticed. For example, we observe that OPLS-AA systematically overestimates the 
energy of 3io-helical structures in vacuo. To exclude any unwanted impact of the overestimation 
of 3io-helices, we performed additional constrained basin-hopping searches in which we forced 
certain H-bonds of the molecule to remain 3io-helical, again followed by individual, unconstrained 
PBE+vdW relaxations. 

The PBE+vdW relaxed conformers were sorted into "families" according to their H-bond pattern. 
We define an H-bond to be present when an O acceptor atom is closer than 2.5 A to a H donor atom. 
Within each H-bond family thus defined, small conformational variations are still possible, e.g., by 
slight bends of the backbone atoms or different rotamers of the LysH+ side chain. Typically, the 
lowest-energy PBE+vdW conformer is found among the family members arising within 5 kcal/mol 
(0.2 eV) of the lowest-energy force field conformer. 

Harmonic vibrational frequencies and intensities were computed from finite differences. The accu- 
racy of the vibrational frequencies calculated in this manner is estimated by analyzing the rotational 
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and translational modes. These modes should be at zero frequency, and thus the deviation ob- 
served gives a hmit for the accuracy of the rest of the frequencies. We have not observed deviations 
of more than 2 cm~^ in any calculations. Harmonic free energies were calculated in the harmonic 
oscillator/rigid body approximation. Vibrational density of states beyond the harmonic approxi- 
mation were calculated from the Fourier transform of the velocity time autocorrelation function, 
taken from ab initio molecular dynamics trajectories. 

For all molecules containing alanine residues and for the Li+ terminated model structures 
discussed in this work, no extensive conformational search was performed. These peptides are 
structure models used specifically for a computer experiment to determine the development of 
low-frequency vibrational modes with increasing helix length for two different terminations. Their 
geometries are fully relaxed PBE+vdW structures. 

For additional details about these calculations we refer the reader to the Supporting Information. 

III. RESULTS AND DISCUSSION 
A. Conformational energy hierarchy 

Fig. [l] summarizes the energetic ordering of the lowest-energy (PBE+vdW) H-bond families for 
n=4-8, found employing our search strategy. Only the energy of the lowest energy structure 
belonging to each family is reported, and families are included up to 0.12 eV 3 kcal/mol) of the 
lowest identified minimum of the PBE+vdW PES. a-helical conformers are highlighted in red. The 
3io-helical conformer for n=4 is highlighted in blue. We define purely a- (or purely 3io-) helical 
conformers as those where, counting from the N-terminus, all the backbone CO groups at residues 
i are either connected to NH groups at residues z + 4 (or z + 3) or to the LysH+ side chain (usually 
the final three or four CO groups at the C terminus). Coordinates and a more detailed analysis of 
all the geometries shown in Fig. [l] are given in the Supporting Information. 

The conformer associated with the lowest-energy PES minimum for n=4 is rather small, connecting 
almost all its backbone CO groups to the LysH+ termination. The remaining H-bond at the N- 
terminus is bifurcated, comprising an a- and a 3io-helical H-bond. This conformer could therefore 
be classified as the smallest possible a-helical prototype in this series. In contrast, the lowest- 
energy PES minima for n=5 and 6 are not simple helices. Each contains an "inverted" H-bond 
where one CO group points to the N-terminus and its connecting NH group points to the C- 
terminus, producing more compact structures. In Fig. [TJ these bonds are highlighted in yellow 
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FIG. 1: Energy hierarchies (thick horizontal bars), obtained with the PBE+vdW functional, for the 
low-energy H-bond families of Ac-Alan-LysH+, n=4-8. For each n, we include the conformer 
representatives of the lowest-energy families up to 0.12 eV from the global minimum, as defined by local 
minima of the PBE+vdW potential-energy surface (PES). We also show their hierarchy after adding the 
corrections for the harmonic vibrational free energy Fharm. at T=300 K. Numbers labeling each family, as 
well as colored structural representations are shown. The placement of the structure pictures in the shaded 
gray areas is not directly related to the energy axis {y axis), which applies strictly only to the horizontal 
bars. For these, the a-helical conformer of n >5 is chosen as the reference (zero) energy. For n=4. Family 
1 is taken as the reference, a-helical conformers are highlighted by red bars. The 3io-helical conformer for 
n=4 is highlighted by blue bars. The shaded areas in the energy hierarchies indicate the energy difference 
between the a-helix and the nearest non- helical conformers. 

and pointed to by an arrow. For n=5, we have previously denoted this conformer as "g-1" [46j. 
For n=7 and 8, the lowest-energy PES minima correspond to a-helices. In each case, they are 
closely followed by a conformer that we characterize as compact /globular (Families 2 of n=7 and 
8, with an energy separation of 20 and 60 meV respectively). Thus, we already observe a cross-over 
with peptide length to a-helical lowest-energy minima of the PES at n=7. However, based on the 
energy hierarchies of the structures of the local PES minima alone, one would not expect a purely 
helical ensemble of conformers at room temperature at n=7 or 8. Simple Boltzmann factors would 
indicate a mix of structure candidates. Yet, the experimental work by Kohtani and Jarrold [1^ 
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does indicate a complete room-temperature structural cross-over at n=8 at the latest, albeit based 
on a completely different line of reasoning (water adsorption behavior of "helical" versus "globular" 
conformers) . 

For the low-energy conformers in Fig. [T] we also compute and show in the same figure the impact 
of the vibrational free energy at room temperature (T=300 K) in the harmonic approximation [74j . 
Remarkably, the relative stability of the a-helices is systematically enhanced by the vibrational 
free energy contribution for all n. For n=5 and 6 the a-helical conformers move down in energy 
with respect to the (non-helical) lowest energy PES minimum. For n=7 and 8 the a-helices now 
become the isolated minima. In detail, we observe: 

- For n=8, the energetic interval between the a-helical lowest energy conformer and the nearest 
globular one (red shaded areas) now amounts to 0.14 eV. The additional Family la at 0.13 eV is 
another a-helix with a slightly modified terminating H-bond network. With the vibrational free 
energy included, the energy hierarchy is thus consistent with the experimental claim that a-helices 
dominate over all other possible conformers in gas-phase experiments for n=8 [15, 48j. 

- For n=7, the same qualitative picture emerges. Here, the next remaining conformer (Family 3) 
at 300 K is a mixed 3io/«-helix. The competing compact conformers (representatives of Family 2 
and 2a) are significantly destabilized by Ff^arm.iT)- 

- For n=6, the a-helix emerges as the room-temperature minimum free energy conformer, but the 
competing non-helical PES minimum. Family 1, remains close (50 meV). 

- For n=5, the difference between the a-helix and Family 1 (g-1) decreases by 60 meV at 300 K, 
compared to the PES minimum, but Family 1 remains overall more favorable. 

In essence, there is an uniform stabilization of helical over more compact (globular) structures as 
the temperature increases. The stabilization tendency increases with peptide length, confirming 
quantitatively and systematically from first principles the related observations in Refs. |T3 | [52 t [53] . 
We next demonstrate that it is indeed the vibrational entropy term which is critical, and then 
pinpoint the physical reason among the low-frequency vibrational modes. 

B. Origin of the entropic stabilization 

In Fig. [2] we analyze the individual quantities composing the vibrational free energy differences 
between: Families 1 and 2 of n=4. Families 1 and the a-helical conformers of n=5 and 6, as well 
as Families 1 (a-helical) and 2 for n—7 and 8. The energy terms plotted are the PBE+vdW PES 
minimum energy, the harmonic internal energy (containing the zero-point energy) AUharmX'^)^ 
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and the entropy term TSharmX"^)- They are reported in Fig. [2] as energy differences, taking 
the non-hehcal conformer of each n as the reference, such that negative slopes correspond to a 
stabihzation of the a-hehces. Upon inspection of Fig. [2j we observe a monotonic stabihzation of 
ah hehcal conformers with increasing T. While for the shortest molecule (n=4) there is hardly any 
observable effect, the stabilization trend is enhanced with increasing length. For n=6, we predict 
a cross-over of the lowest-energy structures at T ^ 150 K. It is clear that among the individual 
contributions to the vibrational part of the free energy Fharm{T), the entropy term TSharm{T) is 
indeed always the most important helix-favoring term. The zero-point-energy (AUharm r=0), 
is also favorable, but on a smaller scale. 
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FIG. 2: Energy differences for each term of the harmonic free energies, as a function of temperature, for: 
Families 1 and 2 of n=4; Families 1 (g-1) and 2 {a) of n=5; Families 1 ("g-1" like) and 3 (a) of n=6; 
Families 1 {a) and 2 (compact) of n=7; and Families 1 (a) and 2 (compact) of n=8. In each case, the 
non-a- helical conformer was taken as the reference (i.e. negative slopes mean a- helix stabilization). The 
PBE+vdW energy differences are plotted with a yellow dash-dash-dot line, the harmonic internal energy 
AUharm. as black dash-dot lines, the entropy term TSharm. as red dashed lines, and the sum Fharm. as 

green full lines. 



The observed stabilization of the a-helical conformers can be understood by analyzing the lowest- 
energy vibrational modes of different conformations. For simple helices, these modes have been 
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TABLE I: Position of the lowest vibrational mode, in cm~^, for the lowest energy conformers of each 
family shown in Fig. [T]for n=4-8. We use red characters to indicate the a-helical conformers, a * symbol 
for the g-l-like conformers of n=5 and 6, and blue characters for the 3io-helical conformer. The numerical 
accuracy of all computed frequencies is 2 cm~^ or better (See supplementary material). 
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estimated and analyzed before [65j. The important point here is that we have available a direct 
comparison against many other, non-helical structure candidates. These modes dominate in the 
harmonic free energy expression at room T, so that conformers with lower-frequency modes will 
effectively be stabilized over others. In Table [l| the frequencies (in cm~^) corresponding to the first 
normal modes of all conformers shown in Fig. [l]are reported. The lowest-energy a-helices, marked 
in red in Table [ij show first vibrational normal modes between 8 cm~^ and 13 cm~^. The same 
is true for families la of n—7 and 8, which only differ in details of the termination. In contrast, 
the conformers that have first vibrational modes of at least 20 cm~^ are all compact non-helical, 
for example: 28 cm~^ for Family 2 of n=7, and 22 cm~^ for Family 1 (g-1 motif) of n=5. The 
g-1 motif has the first vibrational mode lying around 20 cm~^ for n=5, 6, and 7 (marked with a * 
symbol in Table [l]). 

Of course, the energy destabilization of the compact conformers with respect to the helices will 
depend not only on the first vibrational mode but also on the overall distribution of low-frequency 
modes. To illustrate this point, we show the frequencies of vibration lying between and 50 cm~^ 
for Families 1 {a) and 2 (globular) of n—1 and 8 in Fig. |3]^a) and (b). The a-helices have lower 
first vibrational normal modes, and have slightly higher (one or two modes) densities of modes in 
this region than the globular conformers. 

Here, a comment regarding structures that are more elongated than a-helices is in order. It is 
customary to compare the stability of a-helices with 3io helices and /3-sheets or fully extended 
structures (e.g., in references [25 ] [26 1 1661472] and many others). Following our rationale above, the 
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more extended the structure is, the softer the low vibrational modes will be. This indeed happens 
for the most extreme case, the fully extended structure (FES), where we find the first vibrational 
mode to lie around only 2 cm~^ (calculated for n=8 and 15). In fact, entropically stabilized /3- 
sheets in neutral polyalanine in the gas-phase have been suggested in Ref. [33j. For 3io helices, 
in our own structure searches for n=4-8, we find the first vibrational mode to lie always very 
close to their a-helical counterpart. Accordingly, it is the enthalpic energy difference that favors 
a-helices specifically over 3io for all molecules studied. For infinite periodic structures, calculation 
of phonons and vibrational free energies for a-, 3io-, and vr-helices, and the FES in Ref. ^6] 
corroborate our results. There, all helices are destabilized with respect to the FES at 300 K. The 
TT-helix, which is the most compact among the helices studied in Ref. [66j, is most destabilized by 
the vibrational entropy term. 
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(b) n=8, harmonic 
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(c) n=8, anharmonic, 300 K 
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FIG. 3: Vibrational normal modes with frequencies up to 50 cm~^ for (a) Family 1 (a) and Family 2 

(globular) of Ac-Ala7-LysH+ and (b) Family 1 {a) and Family 2 (globular) of Ac-Ala8-LysH+. (c) 
Vibrational density of states taken from the velocity autocorrelation function of 21ps of microcanonical 
AIMD at (T) = 300 K, for Family 1 (a) and Family 2 (globular) of Ac-Ala8-LysH+. Light colored bars in 
(c) correspond to the harmonic normal frequencies of vibration [same as (b)]. 
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1. Beyond the harmonic approximation: Ah initio molecular dynamics 

We next show that the observations discussed in the last paragraph should hold also in a fully 
anharmonic picture. Even at relatively low temperature, where the structure still stays generally 
constant, we expect first the inherent anharmonicity of the local, nearly harmonic PES, then the 
lowest-barrier transitions between basins (side chain rotations), and then transitions between lo- 
cally different backbone conformations and H-bond networks [El [73] to contribute. Unfortunately, 
a direct calculation of these terms (e.g., by thermodynamic integration) is computationally pro- 
hibitive in DPT. We can, however, use explicit ah initio molecular dynamics simulations to gain 
some msig hts. In Fig. |3j;c), we compare the Pourier-transformed velocity autocorrelation functions 
of the Pamily 1 (helical) and Pamily 2 (compact non-helical) conformers of n=8, extracted from 
explicit microcanonical ah initio molecular dynamics simulations (21 ps total time, 1 fs time step, 
initially thermalized to approximately room temperature). The Pourier transform of the velocity 
autocorrelation function corresponds to a vibrational density of states (VDOS). At T—{) and for 
classical nuclei, the autocorrelation function should reflect only the harmonic eigenmodes of Pig. 
[Sj^b). Compared to these modes, the onset of the calculated VDOS at T^SOO K is noticeably 
shifted towards lower frequencies for both conformers, but the onset frequency for the helix is still 
significantly lower than for the non-helical structure. Thus, the lower vibrational frequencies of 
the helix are also carried over to the full (anharmonic) motion of the conformers. In addition, the 
integral over the VDOS up to 50 cm~^ is 4% larger for the a-helical Pamily 1 than for the compact 
Pamily 2, i.e., the general downshift of frequencies in this region is also preserved. 
Regarding the overall local structure stability during the AIMD simulation, we observe that the 
H-bond pattern of Pamily 2 (compact non-helical) stays essentially the same throughout the entire 
simulation (see Pigure S.l in the Supporting Information). In contrast, Pamily 1 displays local 
structural fluctuations in the helical part, occasionally forming short-lived 3io-like H-bond con- 
nections. Similar fluctuations also occur in simulations of longer helices (e.g., Ac-Alai5-LysH+ in 
Ref. [36j). It seems plausible that the overall greater "floppiness" of the helix compared to the 
more compact non-helical H-bond network adds another favorable entropic contribution at room 
temperature. This direct observation that we make has been guessed as the motive for the loss in 
entropy observed in a-helical formation, compared to a FPU helical structure in Ref. [70j. Ac- 
cording to our reasoning, the FPU structure, being less compact than the a-helices should indeed 
show more fluctuations, just like a-helices do if compared to more compact conformers. 
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C. The role of the termination 

Finally, we show that, especially for the short peptides considered here, there are two independent 
aspects to the stabilization of helical conformers: length and termination. To isolate their role, 
we examine the character of the lowest-frequency modes as a function of peptide length (also for 
longer helices) for two different terminations. The first is the LysH^-terminated series, which is 
the main subject of this work. We contrast this series with Li+-terminated polyalanine molecules 
Ala^i-Li^, a much more rigid termination as we shall see. For the latter, we assume a-helices for 
all n (fully relaxed in DFT-PBE+vdW), which are at least locally stable when Li+ is placed in 
contact with the last three residues at the C terminus. In Fig. jij^a), we compare the position of 
the first and second vibrational modes of these a-helical geometries of Ala^-Li+, n=5-15, and 20, 
with the a-helical geometries of Ac-Ala^-LysH+, n=5-15, and 19. For both peptide series, there 
is a monotonic decrease of the first and second vibrational frequencies for n >7, but the starting 
point is much higher for the Li+ termination (26 cm~^ at n—7) than for the LysH+ termination. 
A softening trend of the respective modes in neutral polyalanine helices with increasing length has 
also been observed in Ref. [65]. In order to characterize this vibrational mode in more detail. 
Fig. (ij^b) visualizes the displacement of the backbone atoms of Ac-Alaig-LysH^ when deforming 
this molecule along the first vibrational mode. Fig. Qc) shows the relative length changes of 
the hydrogen bonds in the structure upon deformation along this mode. Subfigures (d) and (e) 
show the equivalent data for Ala2o-Li^. For both molecules, the vibration spans the helical part of 
the structures. For the LysH+ terminated molecule, we show that the actual LysH+ termination 
(connected to the last four CO residues) is clearly involved in the vibration. In contrast, the 
hypothetical Li+ charged termination constrains especially the C terminus to be much more rigid, 
as evidenced by the almost zero change of all Li-0 distances. Here, the N-terminus is much more 
involved in the lowest-frequency modes. The same trends are observed for the smaller molecules 
in both series. Movies containing 3D visualizations of the first vibrational modes for helical and 
compact structures are contained in the Supporting Information. We thus conclude two points: 

(i) Helices are entropically favored by allowing delocalized, soft low-frequency modes that we do 
not observe in competing, more compact conformers of the same LysH+ termination (evidenced 
by Table [l]). 

(ii) For short conformers, the already electrostatically favorable LysH+ termination [291 132] is 
additionally helpful by allowing soft, delocalized modes to include the termination also for short 
conformers, in contrast to the hypothetical, much harder charged termination by Li+. For long 
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enough helices these softer low frequency modes should exist regardless of the termination. 

IV. CONCLUSIONS 

In summary, we show from first principles, quantitatively, and for a particularly well studied series 
of polyalanine peptides how helices emerge with length and temperature as the leading structural 
pattern from a vast array of possible competing conformers. The crossover to helical stability 
with length is already apparent based on local structural minima of the potential-energy surface 
alone, due to the critical role of H-bond networks including their cooperativity [231427] as well as 
that of vdW terms [27j. In addition, the contribution from softer low- frequency vibrational modes 
acts to stabilize helices with increasing temperature over their more compact competition. The 
specific experimental claim [15j of exclusively helical conformers at and above n=8 is thus explained 
by both enthalpic and entropic effects acting together at finite temperature. Ab initio molecular 
dynamics simulations corresponding to approximately room temperature suggest that these trends 
are further strengthened by anharmonic effects. 

The emergence of room-temperature helix stability with length in Ac-Ala^-LysH+ is thus the result 
of a subtle balance of enthalpic and entropic terms. In a non- vacuum environment, further terms 
would obviously contribute, but we expect the fact that helices, in general, allow locally softer 
vibrational modes to hold. Regarding the overall ubiquity of the helical motif in folded peptides 
and proteins, at least, we here show that low- frequency modes will be a significant quantitative 
contribution. 
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FIG. 4: Characterization of the lower vibrational modes of helices considered in this work, (a) Positions 
of the first and second vibrational modes of a-helical geometries of Ala^-Li^, n=5-15, and 20 (red) and 
a- helical geometries of Ac-Alan-LysH+, n=5-15, and 19 (black). The accuracy of these numbers is of 
around 2 cm~^. (b) Displacement caused on the backbone atoms of Ac-Alaig-LysH^ by deforming this 
molecule in the direction of the first vibrational mode, with red shaded areas corresponding to positive 
displacements and green shaded areas to negative displacements (hydrogen atoms not shown), (c) Relative 

change in H-bond distance when displacing 1 (normalized) unit of the first normal mode for 
Ac-Alaig-LysH+.We number the carbonyl groups from the N-terminus to the C-terminus. (d) Same as (b) 
for Ala2o-Li^. (e) Same as (c) for Ala2o-Li^ (in this case the last three carbonyl groups do not make 

H-bonds, but are connected to the Li+ ion). 
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modes of helices and compact conformers. 
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